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In their Comment, Greer et al (i) put us in charge of a pretended wrong claim, which we never 
made in Phys. Rev. B 78, 115315 (2008), where we criticized a method (DG) proposed by two of 
them, (ii) incorrectly claim that the DG method can reproduce the conductance quantum go, but 
(iii) to deduce go for a toy model, they carry out calculations within the standard Landauer method, 
which has nothing to do with the DG's. We present results for their model obtained within the 
DG method, which demonstrate that the DG method fails as lamentably as in the examples we 
presented in our earlier work. We also analyze the physical reasons why the DG method fails. 

PACS numbers: 72.10Bg, 72.90.+y, 73.23Ad, 73.63.Kv 73.23.Hk 

I. INTRODUCTION 

In their Comment,^ Greer, Delaney, and Fagas (henceforth GDF) claim that they can resolve the conundrum 
expressed in our paper— (referred to as I below), where we challenged a method (DG) put forward by two of them in 
Ref. y and subsequently used in their Refs. |j-|6|. Unfortunately, their attempt to refute our criticism of I (a reference 
which GDF even omit in their bibliography list) is based on a series of inaccuracies and statements without any 
support. This starts from the abstract of Ref. wherein GDF incorrectly state that we criticized a so-called MEGS 
(many-electron correlated scattering) approach, which was not even mentioned in I, nor defined in the works yHy that 
made the object of our criticism. Before discussing the controversial issues in detail, we emphasize the following: 

(i) Using the acronym MEGS, GDF obviously distract attention from the fact that the (DG) method^^^-S envisaged 
by us is manifestly a variational one. In their "brief introduction to the MEGS approach" (Sect. II of Ref. |l|), they 
completely omit to state that the variational ansatz (namely, the minimization of the total energy in the presence of 
an external bias V) is an essential ingredient of the DG method; in that introduction, they do not define a transport 
approach, but merely review well known properties of the Wigner function (WF) f{q,p). 

(ii) With this acronym and "carefully" chosen quotations from I, by only mentioning in passing "the variational 
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structure of the MECS calculations performed to date ..." (Sect. V), without noting that all their calculations 
criticized in I are variational, they aim to convey the false impression that the DG method merely consists of open 
boundary conditions (OBCs) expressed in terms of the WF, and that I only criticizes these WF-OBCs. We stress 
from the very beginning that what we criticized in I is the DG variational approach as such, and in particular the 
WF-OBCs in the specific context of that method employed to a finite cluster; the DG approach, as any transport 
theory, cannot be reduced to whatever boundary conditions (BCs). 

(iii) None of the results of GDF (their Sect. IV) were obtained within the DG method, but simply within the 
well-established Landauer approach for an uncorrelatcd toy model (GDF model hereafter) ii 

GDF mainly claim that: our critique in I of the DG method is invalid, and they present calculations within a simple 
model that reproduce the conductance quantum go- The first claim is wrong, as explained in Sect. |lTl In Sect. Illli 
we show that the second claim is irrelevant, because their derivation of go is done within the Landauer approach 
and not within the variational DG approach criticized by us. Some challenges to the results thus obtained by GDF 
are presented in Sect. IIVI Next, in Sect. |V] we present results for the GDF model deduced within the DG method, 
which clearly demonstrate the failure of that method. The other issues raised by GDF are addressed in Sect. IVII In 
Sect. IVlJ we indicate the physical reasons why the DG method fails. Conclusions are presented in Sect. IVlIIl 

II. REFUTING THE GDF CRITICISM TO OUR WORK 

The main presumptions made by GDF to refute the criticism of I are that: (i) we incorrectly presumed an asymmetric 
injection of momentum, (ii) we incorrectly assumed that there is no chemical potential difference when applying OBCs, 
and (iii) we appear to be confusing in the application of OBCs when using either momentum or energy distributions. 

All these claims are wrong. To arrive at our basic Eqs. (19-21), which yielded the results presented in Figs. 2 — 
5, and 7, demonstrating the lamentable failure of the DG method, we made exactly what that method prescribes: 
energy minimization, Eq. (1), with WF-OBCs, Eqs. (5, 6), wave function normalization, Eq. (9), and, if not otherwise 
specified, current conservation, Eq. (8). More specifically: 

To (i): A momentum asymmetry was assumed neither in the analytical results of Eqs. (19-21), nor in the numerical 
results of Figs. 2 — 5, and 7. In our WF-OBCs, Eqs. (5, 6), the "momenta" pl.r hi the LHS and RHS are the same. 
Otherwise, we should have written, e.g., f{qL.,R,PL,R± Sp) = foiqLM,PL.,R), Sp = 0{V). 

To (ii): This is wrong, we do assume a nonvanishing chemical potential difference (eV). The bias V enters our 
calculations, because the quantity to be minimized of Eq. (1) does depend on it via W. Consequently, the wave 
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function ^ obtained within the DG procedure does contain an asymmetry between the left and right electrodes for 
V ^ 0, and this asymmetry reflects itself in the WF at the boundaries. Without this dependence, e. g., the quantity 
(5/ = / — /o 7^ for electrons flowing from the device into electrodes displayed in Fig. 2 would have been zero. 

To (iii): By faithfully applying the DG method to well-defined models in I, we were simply faced with a mathematical 
problem whose solution is exact and unique within the linear response theory, and unambiguously yields unphysical 
results. We needed not explicitly discuss, e. g., the distinction between momentum and energy distributions. 

The numerical results presented in I, which arc completely unphysical, arc nothing but the results of the variational 
DG method, because in I we did nothing else than exactly what that method prescribes. Therefore, any critique to I 
represents in reality a critique of the DG method itself. 

III. REBUTTAL OF THE CLAIM ON THE REPRODUCTION OF THE CONDUCTANCE QUANTUM 

GDF write: "/n practice, the" WF "is used to constrain the momenta flow out of the electron reservoirs and into 
the scattering region" (Sect. II). This "practice" does not obviously refer to the Comment; simply, GDF use Eqs. (3-5) 
nowhere in calculations (not even in their envisaged works). Relevant for the present debate would have been to apply 
the variational DG approach (including its original WF-OBCs) and see whether it is able to correctly describe the 
transport, and not to show that another (Landauer's) approach can reproduce correct results (conductance or WF). 

In the Comment, GDF do nothing else than apply the standard Landauer approach to deduce the conductance 
[see their Eqs. (12) and (13)] and to compute the WF in electrodes' middle. So, they cannot pretend that the 
variational DG approach correctly reproduces the conductance of their simple uncorrelated model. Their statement 
in Introduction that in "Seci. IV a calculation of conductance quantization . . .is given using the MEGS construction" 
is not true. Equally wrong is their claim that the BCs "as formulated in MEGS applied to" their "model reproduces 
the well-known result of conductance quantization" . We did not challenge the Landauer approach, and the fact that 
using it GDF can deduce a correct result {go = e^/h) for an uncorrelated model, which is trivial from the point 
of view of the Landauer approach, is not at all astonishing: the transmission through a barrier of vanishing height 
{y — >■ 0) is equal to unity. This correct result, deduced within a correct approach, has absolutely no relevance for 
the validity of the variational DG method, and our critique of the latter is not in the least affected. In the GDF 
calculations, BCs are imposed not by using the WF but by simply postulating standard scattering forms of the 
asymptotic single-electron wave functions [Eqs. (7) and (11)]. They use these forms, in a kind of a posteriori check, to 
show that the WFs computed within the Landauer approach, in and out of equilibrium, at carefully chosen boundaries 



behaves as expected physically. Computing the WF within the Landauer approach (and not within the DCs) and 
claiming that ^^indeed it is observed that in the single particle case constraining the incoming momentum inflow via 
the Wigner function implies solving the Schrddinger equation ..." (Sect. IV. A of GDF) is totally unfounded. The 
GDF calculation represents in itself a logical inconsistency hard to understand: within Landauer calculations, they 
assume Fermi distributions (shifted in energy by /i/j — /i^ = eV for V ^ Q) va electrodes, and show that the WFs 
(approximately) just look like the Fermi functions already postulated. Within this philosophy, one can solely infer 
that the Landauer approach is indeed self consistent, but this is trivial and irrelevant for the validity of the DG's. 

As noted in Introduction, GDF wish to reduce the whole debate to the WF-OBCs, but no transport theory can be 
solely reduced to a certain type of BCs. In I we criticized the WF-OBCs in the specific context of the variational DC 
method, and did not challenge the WF-OBCs in general. See the end of the third paragraph on page 2 of I and the 
lines after Eqs. (3) and (4) in I], or what we explicitly wrote (last paragraph of I): "While the idea of formulating 
boundary conditions in terms of the Wigner function for correlated many-body systems is interesting, the manner in 
which it was imposed in" Ref. y "turns out to be inappropriate. " 

By deducing a simple result within the Landauer approach and not within the DG's, GDF raise very serious doubts 
that they can also derive the correct result within the latter. Confirming these doubts, we shall demonstrate in Sect.lVl 
that for the simple model of their choice, the DG approach fails as lamentably as in the cases presented in I. 



IV. CHALLENGES TO THE GDF CALCULATIONS 



The CDF's Landauer calculations affect in no way our critique to the DG method. However, especially because 
GDF attempt to present these results as if they were related to Ref. 3 , we note several errors of Ref. 

For the GDF model, the transmission can be trivially obtained analytically. Tie) — '^kLkji/ {kL + kj^)"^, k^ = 



V^me/h, and kn = \/2rfi{e — eV)/h^ and it yields the following exact expression of the current 
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The above ny coincides with the RHS of CDF's Eq. (16). 

Instead of giving these simple formulas, GDF present some qualitative considerations aiming to show that (i) the 
conductance quantum is correctly reproduced, and (ii) this result has something to do with the results and the electrode 



sizes {L w 2nm) of their controversial works.'^^— To deduce go in Sect. IV. B, they approximate T(e) ~ T{£f) — 1; 
this amounts to implicitly assume the linear response limit [eV <C £f)- In the last paragraph of Sect. IV. A, GDF 
claim that their considerations to deduce go ^^opply well to electrode lengths as small as 1 nm" L/2). They assume 
that np is large and Afc — 2tt/L is small [this should mean that np — kpL/27T ^ 1 and Afc = 2tt/L <^ kp, 
cf. Eq. (7) of GDF]. For L — 2nm and their completely unphysical value kp = 0.12nm^^ (cf. Fig. 3 of Ref. 
one gets np = 0.0038 (!) and Ak — 3.14nm^^ ^ kp{\). In reality, to mimic gold electrodes (Fermi velocity 
vp ~ 1.4 X lO^m/s, ep — mvp/2 ~ 5.6eV), one needs a fci?-value hundred times larger, kp ~ 12nm^^. Even then, 
one gets np = 3.8, and satisfying the above conditions is problematic. In fact, to derive Eq. (18) from Eq. (17) in 
Ref. Ill GDF need not a large np^ but rather a large ny, which should obviously be much smaller than np. If we 
admit that ny = 2 is a number "much" larger than unity and "much" smaller than the other "large" number ^ 4, 
Eq. (P) leads to V ^ 1 volt. 

How poor is the linear approximation in this range and how bad is the description based of electrodes with L = 2 nm, 
exactly the linear sizes of the Auis and Au2q clusters used in GDF's ab initio works, '^"^ can be seen in Fig. \T\ and 
any further comments are superfluous. If, as in the present case, the exact current for L — >■ oo is known, one may 
still argue that the "correct" trend can be unraveled by a "clever" inspection of the curve for very short electrodes 
(L = 2 nm) in Fig. [TJ However, one may more legitimately ask how reliable could be considered a theoretical result 
obtained with very short electrodes in implementations for complex realistic cases (like those of Refs. |3|-l6|) where 
neither the exact values nor the electron spectrum details are known a priori. From the curve for L = 2 nm of Fig. [1] 
one deduces a linear conductance (V — > 0) of ^ lO^^goi i- 6., hundred times smaller than the true value go. 
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FIG. 1: I-V characteristics obtained exactly from Eq. ^ for infinite electrodes with kp = 12nm~^ (as in gold), and for two 
finite sizes L (given in the legend). 



As their results on the WFs obtained within the Landauer approach are irrelevant for the DG validity, we restrict 
ourselves here to a few critical remarks. In Sect. II, GDF mention that the WF, as a function of energy defined in 



the phase space, tends rapidly toward the Fermi distribution with increasing number N of particles in a confining 
potential. Indeed, in their Ref. 9, Cancellari et al showed that J-{ep) ^ ^ _^g^ dq^ dp^ f^^\q^ ,p^) behaves 
as the Fermi distribution f{ep). Notice that the system of that Ref. 9 is confined (while GDF use periodic BCs, 
k„ = inn/ L), and Cancellari et al attribute a physical meaning to a function whose argument is the energy Sp, which 
is unequivocally defined in the uncorrelated case discussed by GDF, and not to a function of momentum p. Contrary 
to them, GDF ascribe a physical meaning to f{qL,R,p), i- e., at fixed locations (without g-integration) and interpret 
p as a physical momentum of electrons, e. g., which move toward right for p > and toward left for p < 0. We do 
not analyze here f{q,p) a.t V OM^ We only point out an error in the Comment for V ~ 0. The exact expression of 
their Eq. (10) deduced from their formulas for finite L is 
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where P = jifpI {fikp)- Inspecting Eq. ([2]), one can immediately see that the RHS vanishes for p = 0. Therefore, their 
formula and the curves for L ~ 2nm and L = 20 nm of their Figs. 3a and b are incompatible (the case qji = +-L/4 
is similar). The Wigner function computed by means of Eq. i. e., GDF's own formula for their sizes is shown in 
Fig. ^u!^ As visible there, the dip around p = [of a width 5p = 2hkF/np oc 1/L, cf. Eq. ([2])], becomes narrow only 
at larger L. If GDF used their formula to compute the WF at ql.r — they would have shown the curve of our 

Fig. [2^, not too much resembling a Fermi distribution for L ~ 2 nm, which mimics their gold electrodes<^ 



2.0 



0.5^ 



(a) 






L=20.26 nm 








\l=2.0S nm 




2.0 



-0.5 -0.4 -0.3 -0.2 



(c) 


li,- \^ 








1=20.26 nm 


















\l=2.0S nm 


-1.5 


1 -0.5 0.5 


1.5 p/p, 



-0.1 q/L 

FIG. 2: Equilibrium Wigner function computed: (a, b) using the GDF formulas, and (c) for a confined system for two sizes L 
(given in the legend) and kp ~ pp/fi ~ 12 nm^^. 



Limitations of the physical meaning of the WF are well knownj^^— Noteworthy, in emphasizing these limitations 
even at equilibrium, textbooks (see, e. g., Ref. |l2|, ch. 3.7, pp. 202-203) choose typical examples just like the one 



invoked by GDF (cf. Sect. IV. A) to claim that the WF is physically useful. The fact that the boundary locations 



whereat the WF was computed by GDF cannot be arbitrary and should be "carefully" chosen (just in electrodes' 
middle to GDF, qL,B. = already raises warning bells on its appropriateness for applying BCs. Generalizing 

Eq. ^ away from qL,B. = ^L/A by restricting the integration r-range to either electrode in the GDF formulas is 
straightforward. To illustrate the sensitivity of fo{qL,p) to the boundary location, we depict in Fig.[2}D the component 



which defies a simple physical interpretation. 



p = 0, a classical textbook's example (e. g, ch. 8.8, p. 327 of Ref. 
Noteworthy, our results of the genuine DG calculations in I are quite unphysical even with this choice (boundaries 
exactly in electrodes' middle): this is just the choice in Fig. 7, this choice also qualitatively changes nothing in Fig. 3. 

V. RESULTS OF THE VARIATIONAL DG METHOD FOR THE GDF MODEL 

We proceed by showing what GDF should have shown, namely the results of the variational DG method for the 
chosen model. To obtain the linear response of the GDF model, the working equations (5-21) of I (which yield an exact 
and unique solution) deduced from the faithful implementation of the variational DG method without any other extra 
assumption (cf. Sect. |Tl|), can be straightforwardly used. Only a minor adaptation is needed in Eqs. (7), (10), and 
(12), namely to consider a continuous space coordinate —L/2 < x < L/2 instead of a discrete lattice, and this will be 
briefly indicated below. An important strong point of the approach of I is that it allows to determine self-consistently 
the limit of its applicability, namely the highest bias V < Vir compatible with the linear response approximation. For 
this, the expansion coefficients An — 0{V) [4* = Aq'^q + J^n^a ^n'^n, = 1 ~ 0(T^^)] should satisfy the condition 
S = J2n l^nP ^1- As a practical numerical criterion we imposed S — 0.01, which typically yielded Vri ~ Ivolt. 
Therefore, all numerical results presented in this section for V are for V = Ivolt. 

Without applied bias (V^ = 0), an eigenstate k of each electron out of the N nonintcracting electrons considered is 
characterized by a wave function 4>k{x) and an energy satisfying the Schrodinger equation 

^2 £ 
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<i3k{x) = ek4>k{x). (3) 



2m d x"^ 

The second quantized electric current i{x) and Fano F{q,p) operators needed to compute the linear response read 

eh 
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The action of an applied voltage V{x) is expressed by 



W = y^Vk,k'clckr,Vk.k' =-e j dx(t>l{x)^k'{x)V{x). (6) 
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For the GDF potential, the antisymmetric form V{x) = V/2 if < x < 0, and V{x) = -V/2 if < a; < L/2 is 

more convenient, because it permits to separate the even and odd (^u) many-body eigenstates, as discussed in 
I. For this, it is necessary to further assume a symmetric potential barrier v{—x) — v{x), which can be also included 
to make the calculations a bit more realistic. This (let us call it statical) barrier should not be confused with the 
barrier V{x) related to the applied bias. For the GDF model, v{x) = 0, and — h^k'^/{2m). Within the DG 
approach, the system is confined within —L/2 < a: < L/2, (j)k{±L/2) = 0, and therefore (jjf. (x) ~ ■\/2/Lcos(fc„x) and 
'^'k ~ \/2/Lsin(fc„a;), where fc„ = ttii/ L. For even eigenstates (superscript 17), n = 1, 3, 5, ... is odd, while for odd 
eigenstates (superscript m), n = 2,4, 6, . . . is even. In the ground state ^Po, the highest occupied single electron level 
has k = kp = ttN/L. GDF's electron orbitals (j)k„{x) = exp(ifc„a;)/\/X satisfy periodic BCs, fc„ = 27rn/L, where n is 
a signed integer {—np < n < np). Without asking whether it makes sense to consider systems which are not confined 
within the original variational DG approach, for completeness and for comparison with the exact results for the GDF 
model (Sect. ITV)) . we carried out calculations for both cases (termed confined and periodic below). 

Because the operators W , j, and F{q,p) are bilinear, all their matrix elements (5'n| . . .c\ck'\'^Q) needed as input 
into the working Eqs. (12) — (21) of I can be obtained exactly. All the exact excited many-electron eigenstates ^„ 
that contribute are particle-hole excitations, |^'„) — > \'i>k.k') = (^ l)'^''''°'c^Cfc'|5'o). Here, the sign factor accounts for 
the ordering adopted to fill the Fermi sea. So, exactly as in I, we can present exact full CI calculations done within 
the DG approach, and the results should be exact if that approach were valid. 

Within its variational scheme, the DG approach determines the steady state 5* at F ^ by constraining the WF of 
incoming electrons at the electrode-device boundaries (at ql.b. = Ti/4, as by GDF) to that of the ground state ^Pq 

fiqL.R,±P > 0) = mF{qL,R,±p > 0)1*) = {'^o\FiqL.R,±P > 0)\^o) = foiqL.B.,±P > 0). (7) 

The constrained minimization of the total energy allows the system {—L/2 < x < L/2) to optimize the WF of outgoing 
electrons and therefore the differences below are allowed to be nonvanishing 

Sf{qL,R,TP > 0) = fiqL,R,TP > 0) - /o(<Zl,TP > 0) = {^f\FiqL,R,TP > 0)|*) - {^„\FiqL,R,TP > 0)|*o) ^ 0. (8) 

Due to this fact, by mathematical construction, the DG steady state breaks the time reversal. In Fig. we present 
results for the WF fQ{qL,p) = foiQRTP) in the ground state ^E'q of the confined system. Comparing panels a and c of 
Fig. [5] one can see that, as expected for any physically relevant quantity, the equilibrium WF saturates at sufficiently 
large sizes L and becomes independent on the BCs, let they be periodic or open. Notice that these results are for 
equilibrium {V = 0) and have nothing to do with the DG method. 



Let us now examine the numerical resultsl^ for the differences of Eq. ([S]) computed within the DG approach (Fig. [3]). 
Being for outgoing electrons that are unconstrained, they are indeed nonvanishing. Because Sf{qL,p < 0)[= Sf{qR = 
—QLt^p) for linear response] are real, we note that both the real and imaginary parts of A^'s of I contribute to 
Sf{qL,R,TP > 0)- The latter contribution, denoted by Sfi, does not vanish, as visible in Fig. [3^, which reveals that 
ImAu's related to the current [cf. Eqs. (18) and (21) of I] do not vanish. This holds both for the confined and for 
the periodic case. So, in principle, the DG approach can allow a current flow. The crucial point is now whether the 
electric current associated with the time reversal driven (or, "mimicked" , cf. Sect. lVIl)) by these (5/ ^ is appropriately 
described within the DG approach. Our exact calculations demonstrate that, as in the two cases presented in I, the 
DG approach is invalid: although df{qL.R, TP > 0) 7^ and IniA,; 7^ 0, both in the confined and the periodic cases, 
imposing current conservation or not, the linear conductance gjja computed within the DG scheme vanishes within 
numerical accuracy or, to be on the safe side, gua < lO^^^o for the investigated sizes 2nm< L < 400nm.— 
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FIG. 3: (a) Deviation from equilibrium of the Wigner function at the boundaries Sf{—L/4:,p) = — (5/(+L/4, — p) for the confined 
and periodic cases (superscripts o and p, respectively) [see the main text and Eq. (b, c) The very slow convergence of 

5f{ — L/4:,p) by increasing all the exact multielectronic eigenstates up to high excitation energies Ae — rep {r is given in inset) 
substantially larger than the metallic electrode bandwidth ('^ ef)- Because of rapid oscillations, only the range around the 
Fermi momentum pp is shown in (b). 



In fact, this result is not at all astonishing; on the contrary, it should have been expected. The GDF model is nothing 
but the continuous version of the discrete model of Eq. (22) of I, for the particular choice II — — td^L — td,R, 
El = £r = £g, U — 0. As seen on the SWF-curve in Fig. 3 of I, at resonance {el = sr = £g), the DG-conductance 
also vanishes. In this context, what is worth emphasizing is not that the DG approach yields a(n almost) vanishing 
conductance. In fact, the present calculations confirm the analysis Sect. VI of I that the situation at resonance is 
particularly favorable for the DG approach to predict a vanishing current. But, as shown in I, g^Q can also be 
nonvanishing. Fig. 3 of I depicted situations where g^Q 7^ 0. But, completely unphysically, as seen in that figure. 
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the farther away from resonance, the larger is the conductance goG- Although more tedious, exact DG calculations 
are also straightforward for the continuous-space counterpart of the uncorrelated discrete model of I; to this, one 
should include in Eq. ^ a statical rectangular barrier v{x) — eg9{\x\ — I), wherein 9 is the Heaviside function and 
Q <l < L/Ai^ As far as the exact Landauer approach is concerned, the only difference from the uncorrelated model 
of I is that the Lorentzian decay (cf. Fig. 3 of I) of the exact transmission (Landauer conduction) away from resonance 
there is replaced by an exponential decays with Eg and I here. 

To demonstrate the failure of the DG approach, we have chosen above the same method pursued in I, of faithfully 
applying the DG method. Due to its extreme simplicity, the GDF model also offers another alternative decisive way 
to challenge the DG method. Owing to the fact that the analytical exact scattering solutions of the single particle 
Schrodinger equation are known for arbitrary V , one can pursue a different route. Namely, one should consider small 
deviations of the single electron wave functions with respect to the exact ones {(pk 4>k + and investigate their 
impact on the BG-innciion'^XU = {■^\H+W\^)-u:{^^)-Y.^^^^\,{^\F 

— j'(0))|4'); it should be minimum if the DG method were correct. That is, one should examine 
whether 61-1 = 0{6(j)k) or SH = +0{{6(l)k)'^). The expressions thus obtained are somewhat similar to those worked out 
in our approach of the DG linear response,^ but are obviously no longer limited to the linear response. Moreover, 
one can thus directly scrutinize the validity of the variational approach itself, by constraining instead of the WF 
other properties, which are less sensitive to the boundary locations and, highly desirably, have a clearer physical 
meaning. Thus, one can indeed exploit another advantage of an analysis based upon an analytical model ..." that 
avoids issues associated with linear response approximations, perturbation theory, variational methods in a finite 
basis, specific implementations of electronic structure, or other numerical approximations, thereby allowing a clear 
focus on the physical assumptions made when using the^^ DG "method^^ (cf. Sect. V of GDF). 

The results obtained as described above deserve a separate analysis, which is beyond the scope of this Reply. What 
is important for the present purpose is the unambiguous demonstration that the variational DG method (without any 
other assumptions) lamentably fails even for the model chosen by GDF themselves. 

VI. FURTHER ISSUES 

In Sect. V, GDF claim that we criticized "i/ie use of a configuration expansion to describe transport properties^ 
Such a statement cannot be found in I. GDF should not reduce the transport theories based on a configuration 
expansion (CI) to the DG method. What they actually mean is our clear statement that, even if the DG method 
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were sound, it would be impractical. We showed that even if very many exact eigenstates were included, the current 
within the DG approach would very slowly converge; for a correlated system, the first 300 exact eigenstates out of a 
total of 1225 are insufficient. One may think that the convergence is an issue only for correlated systems. Indeed, in 
correlated systems, for which the DG method was designed, the convergence is extremely poor. 

But let us examine the convergence in the uncorrelated GDF model. Because the current (practically) vanishes, 
let us inspect the quantity > 0) that is directly related to it via time-reversal breaking. To exemplify, 

along with the convergent results, we present in Figs, ^jp and c those obtained by including all exact eigenstates 
with excitation energies Ae below a given value Ae < rep- Both for sizes where ab initio DG computations were 
attempted {L ~ 2nm),'^ and for those [L ~ 20 nm) where they are hopeless, all exact eigenstates with very high 
excitation energies (r ^ 1) much larger even than the metallic electrode bandwidth sp) have to be included to 
reach convergence. Definitely, this were a serious challenge even for trivial uncorrelated systems and even if this 
approach were valid. From the perspective of the severe limitations of the number M of multielectronic configurations 
in nontrivial ab initio calculations, the very slow convergence of the DG-results observed for the GDF model, for 
the two models of I, and many othersiS also becomes an important issue because, e. g., increasing A/" by ~ 50% at 
the limit of feasibility and obtaining a change ^ 3 % in the investigated properties, one can erroneously conclude 
that the results almost converged. Indirectly referring to the poor convergence of the DG method demonstrated 
in I, GDF claimed in Sect. V that ^Hntegrated quantities such as the energy may be better approximated compared 
to local properties such as . . . current density" . But in reliable transport treatments the convergence is needed just 
for the latter. Letting alone that, if current conservation were correctly accounted for, the (position-independent) 
current would also be an "integrated quantity", / = L^^ j dxj{x), let us illustrate the effect of a truncated CI for 
a case relevant just for the electrode size L/2 ~ Inm of Ref. 3. To compare the convergence of the DG method 
to that of a standard calculation, we consider the change in the total energy caused by 7^ in the DG-state, 
AE^^ = (^'|iy|^'), and in the state 5'*^^-' obtained in the first order [©(V)] of the perturbation theory (without any 
WF- and current-constraints) AE^-^^ = {'^^■^^\W\'i'''^^) . By including all exact eigenstsites up to an excitation energy 
Ae = 8eV, we got within obvious notations AE^^^^^ = 0.18A£:^^„ and AsiHrox = OMAE^dlv- The convergence 
is an issue for the DG method, not for a certain particular model. Notice that (i) we considered all aforementioned 
exact eigenstates, (ii) this Ae-value is four times larger than that of the states considered relevant by DG,"^ for which 
DG mentioned (without any detail) an inaccuracy factor ~ 3,^^ and (iii) the convergence dramatically deteriorates 
for a real correlated system. Contrary to DC,— we cannot see any justified manner to evaluate the inaccuracy factor 
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related to a CI truncation from Figs, ^jp and c. 

It is obvious that the variational DG approach does not determine the wave function of an eigenstate (cf. Sect. V 
of GDF), no transport theory should attempt to do this; otherwise, the current would be identically zero unless the 
systems (e. g., superconductors) sustain persistent currents. We can find neither in Ref. y nor elsewhere a mention 
that the central variational DG ansatz would be an approximation. One deduces that the results are exact if one is 
able to consider sufficiently large sizes in full GI calculations. This is the case of the uncorrelated discrete model of 
I. In a valid transport through uncorrelated systems genuinely based on the WF, the current conservation is exact 
(provided that the Fourier completeness is satisfied if spatial and momentum grids are used)»i^ We demonstrated in 
I that the DG method does not automatically satisfy the current conservation, which needs be explicitly imposed. 
We did not misinterpret anything, as GDF argue, we criticized the DG claim that this constraint is necessary only 
for other approaches, which use nonlocal interactions, truncate the molecular orbital basis set or the GI (second 
paragraph of Sect. VI in I), but not for the DG method, as if that method were so good and automatically accounted 
for it. The current conservation is trivial, it is satisfied by construction in the calculations where the current is 
constrained to be position independent, as wc also did in I. It is true "f/iaf . . . considerable care is needed in defining 
finite expansions that are current conserving^^ (cf. Sect. V of GDF, underlined by us), but this does not affect the 
results of I, which demonstrate that without explicit imposition, the DG method violates the current conservation: 
they are deduced within the linear response limit [0(1^)] and full GI calculations. One should not confuse the two 
issues addressed in I: the violation of the current conservation, which is demonstrated by full GI calculations, and the 
very poor convergence, which is demonstrated by studying (as also done above) the effect of progressively increasing 
the number of exact eigenstaies up to values beyond the reach of any feasible ab initio calculations. 

In several places of their Comment, by using the terms "linear response approximation" and "perturbation theory" , 
GDF indirectly mean criticism to I, e. g. in Sect. V, where they mention "f/ie advantage of an analysis based upon 
an analytical model . . . that . . . avoids issues associated with linear response approximations, perturbation theory . . . " . 
GDF should have noted that, as emphasized in Sect. IIVI their derivation of go of Sect. IV. B is also based on the linear 
response approximation: the difference is that they employ it implicitly and heuristically, while our approach uses a 
systematic 0{V) expansion spanning the whole Hilbert space (full GI). 

The current oscillations mentioned by GDF in Sect. V (which are rather irregular fluctuations^'^) have neither to 
do with the simplicity nor with the dimensionality {d — 1) of the models of I. As noted there, for the uncorrelated 
model and the sizes considered by us, the time-dependent density matrix renormalization group (t-DMRG) yields the 
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correct physical result, while the DG method lamentably fails either with or without imposing current conservation. 
Likewise, the correct result is obtained within the standard Keldysh formalismjii despite the fact that the electrodes 
are modeled by the same tight-binding model used in I, and the employed Green functions G{x, x' , e) should have been 
more affected by "nearsightedness" than the density matrix p{x,x'): they depend not only on the difference x — x', 
but also on energy. The exponent d of the density matrix decay given by GDF is wrong; the correct one is {l + d)/2J^ 
GDF should have also noted that, for their Aui^- or ylu20- "electrodes" , each of L/2 ~ Inm, the spatial variation of 
the density matrix within £ ^0.5 nm (Ref. 20) to which they refer is important, and the invoked "nearsightedness" 
problematic. Concerning this, noteworthy, the above >C-value is deduced by calculating p{x,x') in very large systems 
of sizes much longer that £, and not employing short DG-like electrodes, each of L/2 ~ Inm. 

The simple uncorrelated model of I, which is a textbook's example, correctly describes the major features of 
nanotransport, if the well-established approaches mentioned above are applied, but not within the DG method. 
Definitely, the failure is of the DG method itself and has nothing to do with the simplicity of the model. This issue 
regards the GDF model as well, for which their Landauer-type calculations and ours (cf. Sect. II Vp also yield the correct 
conductance, Kohn's principle notwithstanding. Of course, Kohn's principle is relevant for realistic systems, but, 
curiously, GDF did not note that Kohn's tight-binding prediction provides the best overall fit of realistic calculations 
of the density matrix decay, as shown in the work to which they refer (Ref. 20). 

Contrary to the GDF's claim, I is not a comment on the validity of the DG method. In I we simply checked whether 
the DG method is able to describe the most simple uncorrelated and correlated systems. If we wrote a comment, we 
would certainly have raised even more questions than on the fundamental issues of the next section. Letting alone 
minor issues (e. g, the missing value of V in Fig. 1 of Ref. y), we would have asked, e. g., how did DG conceive 
to apply their method to the Kondo effect (cf . last but one paragraph of Ref. 0) . For typical Kondo temperatures 
Tk ^ 1 mK— 1 K, this would imply to handle electrodes longer than the Kondo cloud length, L > ^ fivp / (kBTK) ~ 
10'^ — 10® nni within a CI expansion capable to very accurately describe a huge number of excited eigenstates, especially 
the very low excitations related to coherent spin fiuctuations of energies ^ ksTK ^ 10^* — 10^^ eV. 

VII. DISCUSSION. WHY DOES THE VARIATIONAL DG METHOD FAIL? 



As already stressed in Sect. I, in I we criticized the WF-OBCs in the specific context of the variational DG method 
(which means more than WF-OBCs) and not otherwise. Let us assume that GDF could have demonstrated that the 
WF-OBCs are correct. What would be the implication on the (in)validity of the DG? None, since our calculations 
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were done just by imposing WF-OBCs (i. e., they were assumed to be "correct"), and the fact that the results obtained 
within the DG approach presented in I are unphysical remains unahered. The only implication would be that one 
could formulate more precisely: the DG method fails not because the WF-OBCs per se are wrong, but because, with 
these OBCs, imposing one or more conditions prescribed by the DG to determine the steady state is unphysical. 

Even if (hypothetically) examples could be found where results of a certain approach (in our case, the DG's) were 
acceptable, a single counter-example suffices to demolish it. The results for two examples presented in I as well 
as those for the GDF model of Sect. IVl unambiguously demonstrate the lamentable failure of the DG method both 
for uncorrelated and for uncorrelated transport. This is an irrefutable mathematical demonstration for the simplest 
correlated and uncorrelated, for discrete and continuum-space models, since their derivation uses nothing else than 
the variational DG method prescribes. Consequently, the main objective of the present Reply has been achieved. 

The problem that is really important is not whether the DG approach fails, but rather why it fails. We do not 
present here a comprehensive analysis^iS but for the benefit of a reader interested in this method, we briefiy note the 
following. Besides WF-OBCs, the DG method prescribes the total energy minimization and the usage of a normalized 
many-electron wave function ^ to describe a nonvanishing steady state current. Noteworthy, both conditions refer to 
a finite isolated system. So, within this philosophy, it could be possible to obtain a nonvanishing steady state current, 
i. e., phenomenon that is manifestly irreversible, by merely examining a cluster which is not only finite (and even 

n 

very small, cf. Ref. 3) but also isolated, by imposing certain constraints at certain (very special) locations inside this 
cluster. The DG method uses absolutely no other information than that pertaining to a finite isolated system: this 
system is completely decoupled from the world, and there is absolutely no source of dissipation. At this point, it is 
important to emphasize that, basically, valid approaches to transport fall into two classes: 

(i) Most widely used are transport theories, e. g., based on the Keldysh NEGF or master equations, which also 
consider a finite cluster (wherein possible electron correlations are treated within DFT or more accuratel y^^i^^i^° ), 
but this finite cluster is linked to infinite electrodes. It is this latter ingredient that accounts for irreversibility in 
a physically justifiable manner: the imaginary parts of the embedding (contact) self-energies become nonvanishing 
only in infinite electrodesi^ On the other side, that the DG approach can lead to a nonvanishing current is solely a 
lucky mathematical consequence: certain matrix elements of the Fano operator (not directly related to an observable 
with an unequivocal physical meaning, like e. g., the electronic number operator) computed somewhere inside of a 
finite isolated cluster happen to have nonvanishing imaginary parts, a fact by no means related to a true physical 
dissipation. Therefore, although by no means critically related with the failure of the DG method demonstrated 
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unambiguously mathematically, we reiterate our claim of I, that imposing WF-OBCs within the DG prescription is 
not physically justified. WF-OBCs can be used, but within other approaches (e. g., by solving the Liouville equation 
for the WF - ), for which the present considerations do not in the least apply. Consequence of an ad hoc mathematical 
constraint without a precise physical meaning, the predicted DG current (vanishing or not) exhibits, not concidentally, 
completely unphysical trends (cf. Sect. |Vl and Figs. 3 and 7 of I). Summarizing, in this paragraph we have indicated 
one serious flaw of the DG method. 

(ii) Another class of transport treatments deduces the steady-state current by examining the long time (t — > oo) 
behavior of the wave function at zero temperature (e. g., the already mentioned t-DMRG) or the statistical 
(density matrix) operator p{t) at finite temperatures.^^ In the former, starting from x ^ — L/2 (in the present 
notations) the many-body wave packet is monitored a sufficiently long time t, but before the packet reaches the other 
end {x ~ -\-L/2), since in the absence of any dissipation it will be refiected, a reversed current will appear, and current 
oscillations will last forever. Mathematically, this amounts to derive steady-state properties ( e. g ., electric current) by 



221 explicitly emphasizes 



approaching the limit i — > cx) first, and only then t — )■ oo. The clear physical analysis of Ref. 
the importance of the correct order of these two limits for reaching the steady state; see Eq. (6) there. This represents 
the counterpart of the fact well known in solid state physics on the calculation of the dc-conductivity adc from the 
frequency- and wavevector-dependent conductivity a{q,uj). To obtain the correct result, it is mandatory to take the 
limits q — > and — > (the counterparts of L — ?► oo and t — >■ oo, respectively) in that order. 



adc = lim lim cr(q, w); 

CJ— >0 q— ^0 



see, e. g., chapter 3.8 of Ref. 



121 . The above considerations can be rephrased perhaps in a more direct manner as 
follows. The current operator j does not commute with the hermitian Hamiltonian iJ of a non-dissipative system, 
and consequently 

^^^^[H{t),Jit)]^0. 
A steady state limt_>oo ^(i) characterized by a stationary current 

lim dl{t)/dt = 0, and lim I{t) = I =^ 

t— ^oo i— foo 

can only be obtained because, after averaging the commutator (a nonvanishing operator), ('I'|[H, j]|^')t, the average 
can vanish if one first takes the infinite "volume" limit L ^ oo and then i — > oo. 

In examining a steady state transport, two general principles of the thermodynamics of irreversible processes are 
mostly discussed in the literature: the minimum entropy production^'^ and the entropy maximization.'^* Prigogine's 
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discussion of the steady state within the former principle^S clearly revealed that the aforementioned order of the two 
limits is essential. Within the same principle, the same fact was nicely illustrated in the particular case of the flow 
through a capillary tube connecting two containers with ideal gas at different pressuresi^ Were the system finite, the 
flow would be no more irreversible: after some time, the gas would flow back to higher pressure. In Ref. y, DG claimed 
that they can deduce the variational ansatz, which they used to compute the steady state wave function 5*, from the 



lef.y, 



maximum entropy principle. Except for citing Ref. I26l neither in Ref. 3, nor in later works, ■ or in the Comment any 
detail on this derivation was provided. In Ref. |26| . Jaynes presented quantitative considerations emerging from that 
principle applied to systems of very large number of degrees of freedom based on Shannon's entropy for statistical 
equilibrium; a steady state is not even mentioned. In a subsequent paper cited in a later work by DG,^ Jaynes further 
analyzed the time-dependent case in detail from the perspective of a single basic principle (entropy maximization) 
applied to all cases, equilibrium or otherwise. Again, the case of a steady state was not explicitly considered, and 
consequently the issue of the correct limit order to reach a steady state within the principle of entropy maximization 
was not addressed. Prior to the DG work, there has been attempted to recover the standard Landauer results from the 
maximum entropy inference;^ and it turned out that schemes based on that principle encounter notable difficulties 
even if the limits L ^ oo and t ^ oo are taken in the correct order. Ignoring these difficulties and without validating 
their variational scheme against any well-established theoretical result (as done by us in I), DG put forward an 
approach, wherein, implicitly, the limit order is exactly opposite to the correct one: they consider a time-independent 
wave function ^ (amounting to take t — > oo), and then (by taking the largest cluster it can handle, which is in fact 
very small) attempt to mimic the limit L od. This was the physical reason why in I we employed the acronym SWF 
(stationary Wigner function) for what we now called the DG method. Therefore, within this physical context we do 
not think, contrary to GDF (cf. their Ref. 4), that the term SWF is misleading. Summarizing, in this paragraph we 
have indicated another serious flaw of the DG method, which represents the fundamental physical reason why it fails. 



VIII. CONCLUSION 

GDF seem to have realized that their claim is wrong; in Sect. Ill, they write that our conclusion ^Hhat an asymmetric 
injection of electrons is needed to obtain a current is incorrect, if injection refers to incoming electron momentum 
distributions. . . " (underlined by us). However, the above conditional clause does not apply and therefore our critique 
of I is not in the least affected. 

To conclude, in this Reply we have demonstrated that in their Comment GDF could neither show that the critique 
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is incorrect, nor give even a single example where the DG method can correctly describe a transport property. We used 
their model to complete the evidence on the failure of the DG method presented in I. Letting alone the fundamental 
reasons against this method (cf. Sect. IVII[) . because the failure of the DG method^^— which is a unequivocal math- 
ematical prescription, comprises the simplest uncorrelated and correlated systems described within discrete (Ref. |2D 
and continuum (Sect. |V| spaces, it would be a utopia to presume that real systems could be correctly described. 

Most importantly for readers interesting in the DG method, we have presented not only further exact results showing 
that this method fails, but also indicated the basic physical reasons why it fails. 

DG argued that their essential ingredient, the variational ansatz, is deduced from the maximum entropy principle. 



In Ref. [2^ cited by DGjiS, the author asserted that "if it can be shown that the class of phenomena predictable by 
maximum- entropy inference differs in any way from the class of experimentally reproducible phenomena, that fact 
would demonstrate the existence of new laws of physics, not presently known!"^ This statement may apply to the 
results deduced within the DG method: neither the prediction of the present Sect. |V] (a vanishing conductance on 
resonance), nor a conductance increasing if one moves away from resonance (Fig. 3 of I), or conductance maxima 
of the Coulomb blockade peaks becoming higher and even exceeding go with decreasing dot-electrode coupling td 
(Fig. 7 of I) have been observed so far. These trends are just opposite to those of the available experiments. . . , so 
should one still await the advent of new physical laws in the sense quoted above from Ref. 24? Until then, one should 
not be too surprised if the DG predictions are also at odds with other existing experiments of molecular electronics. 
Until recently,^ although being unable to demonstrate that their theory is sound, DG could claim that their method^ 
produces current values better agreeing with experiment than the other theoretical estimations. The accurate data 
of the beautiful experiment of Ref. |9| have clearly demonstrated that just the opposite is true. To see this, one can 
simply compare Fig. 2a of Ref. 9, Fig. 2 of Ref. y, and, e. g.. Fig. 4 of Ref. 28 among themselves. For instance, at 
V = 1; 1.5 volt, the currents in /iA are ~ 1.7; 6)^ ~ 0.09; 0.13,^ and ~ 6; 9.7.^^ So, without any special fine tuning 
(e. g., contact geometry), a standard NEGF-DFT approach-- can reasonably describe the experimental data, while 
the DG's cannot. 

As a matter of principle, we end by reiterating that whether (un)luckily results obtained within a certain theoretical 
(DG's or whatsoever) method (dis)agree with experiment is neither the only nor the decisive point to assess its validity: 
a comparison with experiment is meaningful only if its physical basis is sound. 
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